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TECHNICAL NOTE D-69I 


STUDY OF SYSTEMS USING INERTIA WHEELS FOR PRECISE 
ATTITUDE CONTROL OF A SATELLITE 
By John S. White and Q. Marion Hansen 


SUMMARY 


Systems using inertia wheels are evaluated in this report to determine 
their suitability for precise attitude control of a satellite and to select 
superior system configurations. Various possible inertia vheel system con- 
figurations are first discussed in a general manner. Three of these sys- 
tems which appear more promising than the others are analyzed in detail, 
using the Orbiting Astronomical Observatory as an example. The three sys- 
tems differ from each other only by the method of damping, which is pro- 
vided by either a rate gyro, an error-rate network, or a tachometer in 
series with a high-pass filter. 

An analytical investigation which consists of a generalized linear 
analysis, a nonlinear analysis using the switching-time method, and an 
analog computer study shows that all three systems are theoretically 
capable of producing adequate response and also of maintaining the required 
pointing accuracy for the Orbiting Astronomical Observatory of ± 0.1 second 
of arc . Practical considerations and an experimental investigation show, 
however, that the system which uses an error-rate network to provide damp- 
ing is superior to the other two systems. The system which uses a rate 
gyro is shown to be inferior because the threshold level causes a signifi- 
cant amount of limit-cycle operation, and the system which uses a tachom- 
eter with a filter is shown to be inferior because a device with the 
required dynamic range of operation does not appear to be available. The 
experimental laboratory apparatus used to investigate the dynamic perform- 
ance of the systems is described, and experimental results are included 
to show that under laboratory conditions with relatively large extraneous 
disturbances, a dynamic tracking error of less than ±0.5 second of arc was 
obtained . 


INTRODUCTION 


For many satellites some form of attitude control is required and, 
in some cases, the control must be fairly precise. The use of inertia 
wheels for precise control of satellite attitude is one of the more prom- 
ising methods within the capabilities of present production and is the 
subject of this report. Various possible methods of damping will be 
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considered, and the destabilizing effect of an integrator, which is used 
to keep the error small, will he examined. After a brief general discus- 
sion, the three most promising systems will be analyzed in detail by means 
of linear and nonlinear analyses, analog computer studies, and laboratory 
experiments. The three systems differ only by the method of damping, 
which is provided by either a rate gyro, an error-rate network, or a 
tachometer in series with a high-pass filter. 

The requirements placed on satellite attitude control systems vary 
from one satellite to another; however, it can be stated generally that 
a reference, or tracking line, in the satellite is required to be main- 
tained in alinement with a specified external reference, or line of sight, 
to a specified accuracy, in the presence of disturbing torques and of 
motion of the external reference . For a stellar or solar pointing satel- 
lite, motion of the external reference will be oscillatory at the orbital 
period, as a result of parallax or velocity aberration. For an earth- 
pointing satellite the external reference will have a constant velocity. 
The most significant disturbing torques acting on an earth satellite will 
probably come from the earth's gravity gradient, the sun's radiation 
pressure, the earth's atmospheric density, and the earth's magnetic field. 

An earth satellite which will suffice as an example for evaluating 
control systems in this report is the Orbiting Astronomical Observatory. 
The most stringent fine control requirements presently contemplated for 
the 0A0 demand that the pointing error be maintained at less than 0.1 
second of arc for one orbital period of approximately 100 minutes. Fur- 
thermore, roll motion about the line of sight must be maintained at less 
than 1 second of arc per second. For stellar sighting the line-of-sight 
velocity will have a maximum of 0.005 second of arc per second, predomi- 
nantly from velocity aberration effects. It is expected that the disturb- 
ance torques will be on the order of 100 dyne centimeters. A brief 
description of velocity and torque inputs is included in appendix A, and 
a more complete discussion may be found in reference 1. 


NOTATION 

Torque constants (dyne cm/ ( radian/ sec) ) : 
B e back emf 

Bp motor friction 

B m magnetic field 

Inertias (gram cm 2 ): 

Jq, satellite body 

Jw 


motor and wheel 


<4- HCO 
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Gains : 

K a motor amplifier, volts/volt 

Ke error amplifier, volts/volt 

Kg rate gyro, volts/ (radian/ see) 

Ki error-rate network, volts/ (volt/sec) 

K m motor, (radians/ sec) /volt 

K s sensor, volts/radian 

Kt tachometer with lead network, volts/ (radian/ sec 2 ) 


outer loop, K s KeK a K m —■ -, (radians/sec) /radian 
gyro loop, KgK a K m y — , radians/radian 


error-network loop, K s KpK a K m -y, radians/radian 


K-t' tachometer loop, Kq-KgK^, radians/ (radian/sec) 

Time intervals ( sec) : 
t time 


t c control time 

t s settling time to within specified error limit 

t sw switching or reversal time of applied motor voltage 

t^cpg incremental error time, or time for error to change from one 

point of zero velocity to the next 

Torques (dyne cm) : 

Tq externally applied disturbance to satellite tody 

Tq_v externally applied disturbance to motor and wheel 

T m motor and wheel output 

^ms 


motor output at stall 


h 

Voltages (volts) : 

V applied to motor amplifier 

V m applied to motor 

Fictitious block diagram voltages (radians/sec) : 

V' K a K m ^ V 

J b 

V K m ^ V m 

Generalized parameters: 

£ damping ratio for second-order equation 

T]_ time constant associated with line-of-sight input s, sec 

t 2 time constant associated with torque inputs, sec 

aj n natural frequency, radians/ sec 

Time constants (sec): 

lag in error-rate network 
T m motor and wheel 

lag in the lead network in series with the tachometer 
Angular positions (radians) : 
cp^ satellite body position 

cp c command or input 

cp e instantaneous error 

cp eQ initial value of cp e 

cp ei linear range of sensor, or error at which sensor output becomes 
limited 

cp ess steady-state error 

Acp 0 incremental error change during t^cp e 


<; j H co 


<; H CO 
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Angular velocities (radians/sec) : 

01 instantaneous satellite tody relative to inertial space 

w r roll about line of sight relative to inertial space 

ojy. instantaneous motor and wheel relative to inertial space 

instantaneous motor and wheel relative to satellite body 
oj^-q initial value of tUy. 

oj 1s line of sight relative to inertial space 

SYSTEM CONFIGURATIONS AND ERROR CONSIDERATIONS 
System Configurations 


A qualitative discussion of system configurations is warranted in 
order to select systems for further investigation. A rather generalized 
block diagram of a control system using a motor and inertia wheel for 
momentum storage is given in sketch (a) . 



General system diagram 


Sketch (a) 

The basic system consists of a sensor, appropriate amplifiers, a 
compensation device, and a motor-wheel combination. An external disturbing 
torque, Tq, is shown acting on the satellite. When a voltage is applied to 
the motor, a torque acts directly on the wheel and reacts directly on the 
satellite body. For simplicity in the remainder of the discussion, the 
motor is assumed to have a linear torque-speed relationship; that is, the 
voltage required for constant torque output increases linearly with speed. 
Such a torque-speed relationship is characteristic of a servomotor. 
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If no compensation is used, the transient response of the system 
will he oscillatory, with the amount of damping determined hy the time 
constant of the motor-wheel combination. For the satellite control sys- 
tems considered, the time constant will he large (since there will he a 
large wheel on a small motor) , so that the system damping will he extremely 
small. Accordingly, some sort of rate device must he used to damp the 
system. Such damping can come from a rate gyro, an error- rate network, 
or a tachometer. The steady-state errors that result from a disturbance 
torque or a line-of-sight velocity may he analyzed hy considering the type 
of compensation used. 

First, consider the effect of a constant disturbance torque, Tq, on 
an uncompensated system. This torque must he canceled hy an equal and 
opposite constant torque, T m , from the servomotor if the satellite is to 
he maintained at a fixed orientation. For an uncompensated system the 
steady torque from the motor will accelerate the wheel at a constant rate, 
and the motor voltage will rise linearly as the wheel velocity increases. 
Thus, in a properly controlled system, a steady-state disturbing torque 
will cause a constantly increasing motor drive voltage and, also, a 
constantly increasing error. 

If, now, damping is added hy means of a rate gyro, the error will 
he essentially the same as for the uncompensated system because the satel- 
lite orientation is nearly constant in the steady state and the gyro signal 
is approximately zero. The same is true for a damping signal from an 
error-rate network. If, however, a tachometer is used for damping, the 
tachometer voltage will increase with wheel velocity. The error signal 
must overcome this additional voltage to maintain control so that there 
is a further, and possibly very large, increase in error. 

Next, consider the effect of a step change in the line-of-sight 
velocity, gj 1s , on an uncompensated system. For proper tracking, the 
steady-state satellite angular velocity, c^, must equal cj 1s . The torque 
to produce the required comes from the reaction of the satellite to 

a corresponding change in the wheel speed, c^. Since, in the steady 
state, the motor voltage increases with Wy., the result of a step input 
of angular velocity of the line of sight will he a change in the steady- 
state motor voltage, V m , and therefore a corresponding change in the 
error, cp e . 

The addition of an error-rate network has no effect on the error 
resulting from a line-of-sight velocity input, since the rate of change 
of the error is zero in the steady state. However, with the addition of 
either a rate gyro or a tachometer for damping, a step change in coi S 
causes a change in gyro or tachometer voltage which adds, perhaps 
significantly, to the error. 

Summarizing, it is noted that damping hy an error-rate network does 
not increase the steady-state errors for either torque disturbances or 
line-of-sight velocity inputs, whereas the rate gyro adds to the steady- 
state error for line-of-sight velocity inputs, and tachometer damping 


CO H -fr- > 
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adds to the errors for "both torque disturbances and line-of-sight velocity- 
inputs. In the case of a stellar-pointing satellite, however, line-of- 
sight velocity inputs are likely to be small, so that the gyro system' s 
deficiency in this respect may not be important. For a system using a 
tachometer, the increase in steady-state error with wheel velocity can be 
compensated for by putting a high-pass network, or lead network, after the 
tachometer in order to pass high frequencies to provide stability during 
transient conditions and to cut off the lower frequencies to reduce the 
error during steady-state conditions. 

For all three systems discussed thus far, the steady-state error is 
constant for a constant line-of-sight velocity input, and is steadily 
increasing for a constant disturbance torque input . The steadily increas- 
ing error looks like a serious handicap, but actually, it is not. As was 
previously shown, this error is necessary in order to force the motor to 
increase in speed to absorb the momentum from the external torque, and the 
error is proportional to the momentum. The motor-wheel unit has a maximum 
amount of momentum which it can absorb; therefore, there is a maximum 
error that can build up . To maintain continuous control of the satellite 
in the presence of a steady torque, provisions must be made for removing 
moment-urn which has been stored in the wheel. (Jet control would be one 
way of doing this.) The system must, then, be designed so that it can 
drive the wheel to maximum momentum without the error exceeding tolerable 
limits. This can be done in the systems previously discussed either by 
having sufficient gain or by adding an integrator between the sensor and 
the motor. 

If the integration is provided in the form of an integrating gyro, 
the maximum torque required from the motor will be much larger than would 
be required if no integrator were present, as shown in appendix B. If the 
integrator is electronic, it will require either a larger motor or a con- 
siderably larger amount of rate feedback, depending upon its location. 

One further possibility for integration would be to use a sampling inte- 
grator which would sample the error at periods which were long compared 
to the natural period of the system, but short compared to the time between 
momentum dumps. Such a system would have essentially no effect on the 
design, except to reduce the forward loop gain required. A reduction of 
the forward loop gain, however, would have the perhaps undesirable effect 
of lengthening the response time. In any case, since the error can be kept 
sufficiently small without integration in any form, it seems reasonable 
to select only the three simpler systems for study. Accordingly, the 
present report is devoted to a theoretical and experimental evaluation 
of the three systems previously mentioned which provide damping with no 
integration of any form. These three systems will be referred to in the 
remainder of this report as the rate-gyro, error -network, and tachometer 
systems; the tachometer system includes a high -pass filter. 



Complete Error Equations 


Error equations for the three systems will now Toe derived. In the 
derivations, the following simplified motor-wheel transfer function will 
he used: 


T m _ KmJw s 

V m 1 + T m s 


( 1 ) 


A complete block diagram of the motor-wheel combination is included in 
appendix C with a brief discussion of the simplifications which led to 
the form used here . 

The rate-gyro system is illustrated in sketch (b) . Here, it is 
assumed that the dynamics associated with the gyro may be neglected. 



Rate-gyro system 


Sketch (b) 

The error equation is derived as follows: 

The open-loop equation for the inner loop is 


“b = 


Jw 

Jb 


K a K m ^ V 


1 + T m s 


^d 

J b s 


Substituting the equation. 


V — KgKgCpe - KgWk 


( 2 ) 

( 3 ) 


and solving for gives 


K 


s^-e^a^-m Te “ j^s ( T m s 

Jw 

T m s + 1 + KgK a K m — 


o>b = 


( 4 ) 


00 H Jpr > 
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Using the equation. 


% = scp c + u ls - scp e 


(5) 


and solving for cp e gives 


J. 


w 


! 

(scp c + w ls ) V s + % + 


A 

Jbs 


<P P = 


s + 


'm. 


s 2 + 


_ | - 

■m r m 


w 


s + 


K s KeK a K m 


‘m 


(6) 


or, expressed with loop gains combined into single symbols as defined in 
the Notation section. 


cp e = 





s + 


Ke' 

T m 


(7) 


The error-network system is illustrated in sketch (c) . 



Error-network system 


Sketch (c) 
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The open- loop equation may he written as 


u>b 


K,-s 

M K e +rWT dS 


K.K, 


J w 


a Jv m j-^ ‘Pe 


1 + %s 


J h s 


( 8 ) 


Using equation (5) and solving for cp e as before gives 




scp + 10 + 


IS J-LS 


s^,I 4 ) s+ 1 


T d T. 




d m 


J- 


w 


J w 


w 


s j + 1 4 + 4 ) s 2 + 1 -±=- + 

1 T d T m 


K s KeK a Km ~ ^ K s K e K a K m — 


m 


T d T m 


s + 


t^t, 


dm 


or, expressed with loop gain symbols, 


(9) 


= 


Td 

*Pc + “is + j^ s 


s^ + 


1 1 
+ 

m. 


T d T. 


S + 


1 


T d T m 


S 3 + ( qr~ + qr" ) S 2 + 

T d T m, 


1 K e ' Ki 
+ -=— + 


T d T m 


' m 


s + 


T d T m/ T d T m 


( 10 ) 


The tachometer system is illustrated in sketch (d). 



Sketch (d) 


03 H -F> 
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The open- loop equation for the inner loop is 


Tm 


KaffmJysV 
1 + T m s 


Substituting the equation. 


V = 


K s K e cp e - 


Kt 

1 + T t S 


T, 


m 


and solving for T m gives 

^ _ KsKeKaKmJys(tts + l)cpe 

m " T t T m s2 + ( T t + T m + Wm* 8 + 1 

Using the equation, 


( 11 ) 


( 12 ) 


(13) 


Tm = Jh s (scp c + w is “ scpe) + T d 


m 


and solving for cp e gives 




scp c + w ls + 


_Tc 

Jbs 


„ . , 1 . 1 . K tK a K m 3 . 1 

S^ + ( h — h ■ ) S + 


T t % 


T+T, 


t T m 


T t T mJ 


KsKeK a K m ~ 

„ , 1 1 KtKaJCm\ 0 / 1 . ” Jh 

s 3 +[■=-+■=-+ - — — s 2 + I + 


T t T 


m 


T t T m 


T t T m 


s + 


KsK e K a K m y 1 
J h 


•m 


T t T m 


(15) 


or, expressed with loop-gain symbols, 


Te = 


Td 

+ “is + jgi 


s2+ (l + l + J^ s+ 1 


T t T m T t T m/ T t T m 


s 3 + 


n 

, a + 

Kt ’ ) s 2 +( 

Vt 

T m 

T tW \ 


1 K e '\ , Ke' 
+ -= — ) s + 


(16) 


T t T m T m J T t T m 
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Steady-State Error Equations 


The steady- state errors resulting from cp c , w ls , and Ta for the three 
systems can he readily calculated If the final value theorem Is applied 
to the error equations just derived (eqs. ( 7 ) , (lO) , and (l6)). 

For a step position input of cp c (s) = cp c /s (with w ls = T,i = 0) , the 
final value or steady-state error is seen to he zero for all three systems, 

since 


‘^ess 

<P C 


= lim scp e (s) 
s-O 


= 0 


(17) 


However, for a ramp position or step velocity input, wis(s) = Uis/s 
(with cp c = % = 0) , the three systems have the following steady-state 
errors: 


Rate -gyro system: 


9ess _ Kg' + 1 
u is K e ' 


Error -network and tachometer systems: 


Tess 1 

Wis Ke' 


( 18 ) 


(19) 


For a step input of torque, Td(s) = T$Js (with cp c = wis = 0) , all 
three systems have a continuously increasing error as a function of time, 
since 


~- f ~ SS = lim scp e (s) -*■ 00 (20) 

Td s-*0 

Although the error due to a step of T^ becomes infinite, equations for 
the increasing error as a function of time can he obtained by taking the 
inverse Laplace transform of the error equations and simplifying them for 
large values of time by setting all terms equal to zero which contain time 
as a negative exponential. The process is very straightforward, except 
that the third-order denominators of equations (lO) and (l6) must be 
factored or else reduced by approximation to second order if tables of 
Laplace transforms are to be used. Reduction of the third-order equations 
to second order is accomplished in the next section of this report. For 
the present, the approximate second-order equations are taken from the 
next section, and the only justification given here is that the steady- 
state error for torque input is dominated by the time increasing term. 


<4 Hod 
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which is not changed hy the approximations. Specifically, for 
T<i(s) = Td/s with cp c = (jj 1s = 0 (and cp e = -cp^) , equations (42), (48), 
and ( 5 6 ) become 


Rate -gyro system: 


Error -network system: 


Tachometer system: 


14 


Ve = 


J b sc 


s + 


•m 


<P e = 


= 


a V 
B + -% 

V 

5 + — 
'm 

^ ( 

, iA 

JbS 2 ^ 

8 + 

_p , Ki 1 

, K e' 

^ * *r 

T m 

T m 

Ta < 
J b^ 2 v 

° + Kt') 


s 2 + 


TtK e ' 

Kt* 


s + 


Kt' 


The resulting steady-state error equations found hy taking the inverse 
transforms of the above equations and simplifying for large values of 
time are: 


Rate-gyro system: 


Error-network system: 


Tachometer system: 


Vess 

T 

1 % K e ' 

+ 

Td 

Jh 

K e ’ 


Vess 

Ki’ 

_ 1 Tm " K e ' 

+ t 

Td 

Jh 

K e ' 



Vess 

1 K t ’ - T t + t 

Td 

Jb K e ' 


( 21 ) 


( 22 ) 


(23) 


These steady-state error equations will he used in selecting system 
parameters as outlined in the next section. 



DESIGN CRITERIA 


Primary Considerations 


The first step in the system design is to select tentatively the 
outer loop gain, the wheel angular momentum, and the motor stall torque. 
The outer loop gain and also the rate-gyro loop gain may he established 
by examining the steady- state error requirements. The angular momentum 
and motor stall torque may be determined if the external torque, control 
time, settling time, and maximum steady- state error are considered. 

The first condition to be examined is the minimum outer loop gain, 
K e ', established by steady-state error requirements. Since K e ' must 
be large enough to keep the steady-state error less than or equal to its 
allowable maximum, equation (19) for the error-network and tachometer 
systems may be rewritten in the form of a requirement on K e ' 


K e ' >- 

c — ■ cp 


essmax 


and equation (l8) for the rate-gyro system may be rewritten as 


, _ w is( K g' + !) 
he > — m 


- q> 


essmax 


An additional requirement is placed on the outer loop gain by the 
error equations (21) , (22), and (23) for torque input. The time for the 
motor amplifier to become saturated in the presence of a disturbance 
torque is the control time, t c . For the Orbiting Astronomical Observatory 
example, t c will be assumed to be equal to the orbital period. For such 
large control times, the increasing error component is dominant in equa- 
tions (2l), (22), and (23), yielding the following approximate equation 
for all three systems: 

a lJL (26) 

Td Jb Kg’ 

This equation becomes a restriction on K e ' when changed to the form 


Kg' > 


^b Tessr, 


For the error- network and tachometer systems, both equations (2b) 
and (27) must be tested to determine the minimum value of Kg' . For the 
rate -gyro system, equation (27) is used to determine the minimum value 


<: 1— 1 co 


00 H > 


2X 
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of Ke* , and equation (25) I s used to determine a maximum ratio of 
Kg'/K e ' "by rearranging the equation and simplifying for Kg' » 1 to 
give 


5bL 

K e ' — w ls 

The simplifying assumption that 

K g ' »1 (29) 

is easily justified if it is noted that in equation (7) the damping term 
in the denominator involves l/r m + Kg'/ T m-> and since r m is large as 
discussed previously, most of the damping must come from the rate gyro, 
so that Kg 'Am » 1/Tm* 

Consideration -will now be given to the angular momentum. The maximum 
momentum which must be stored in the wheel is equal to the disturbance 
torque integrated over the required control time. If the average value 
of T^ is used. 



Jw^Wmax - T dtc 

When equation (30) is combined with equation (27) 


( 30 ) 


K e ’ > 




max 


" J b c Pess- 


max 


( 31 ) 


that is, K e ' must be large enough to produce maximum motor speed without 
the error exceeding its permissible maximum, if all of the available 
momentum is to be used. 


A maximum value of K e ' is set by consideration of the desired 
linear region of transient response for a position step input. If K e ' 
is made much larger than necessary, undesirable limit-cycle operation may 
result. 


Motor stall torque will now be considered. The motor torque and 
corresponding settling time may be roughly estimated for large signal 
response if bang-bang operation is assumed, where full torque is always 
applied in one direction or the other. For a servomotor this means that 
T m is infinitely large. The minimum settling time results if full torque 
is applied in one direction until the output position has reached half 
of its final value and then the full torque is reversed in direction as 
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Response for bang -bang operation 
with minimum settling time 


Sketch (e) 


illustrated in sketch (e) (see ref. 2, 
p . 15) • The minimum settling time, 
t Sm i n , resulting from bang-bang opera- 
tion may be derived as follows: 

The known relation, 


% “ 


Tm 

Jb s2 


(32) 


is used for a step input, T m (s) = T m g/s, 
to obtain 


^(t) 


m m 4-2 

-1 ^ms __ 1 ms T ' 

' Jb^ 3 " 2J b 


(33) 


or 


b J b 

v Tms 


(3^) 


But t = t sw = t Simin /2 for cpjj = cp c /2j therefore. 


:s min 


= 2 



(35) 


Equation ( 35 ) indicates the absolute minimum stall torque required to 
produce a specified settling time. In the design of practical systems 
there are two other considerations which must be included to estimate 
stall torque. First, the actual motor time constant will be somewhat less 
than infinity, resulting in a slightly larger torque requirement for the 
same settling time. The following useful equation for T m is derived in 
appendix C : 


^w^w. 


max; 


*m 


T. 


ms 


(36) 


The second, and most significant, practical consideration is that compen- 
sation is generally insufficient to prevent several overshoots, and the 
stall torque required could well he an order of magnitude larger than 
would be required if no overshoots were allowed. 


o>h 
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Generalized Linear Parameters 


A study can now "be made of the dynamics of the linear, or small- 
signal, response in order to obtain equations for the remaining system 
parameters, and to determine the system linear transient response. 

For this linear analysis an equation must be obtained for cp^ 
resulting from cp c , <jJ 1s , and T^. If the relationship 


A 

k 

1 

8 


<P e = ^c “ % + w is/ s 


(37) 


obtained from sketch (b) is substituted into equations ( 7) > (lO) , and 
(l6) , and the resulting equations solved for cp^, the following equations 
are obtained: 


Rate -gyro system: 


( Pb = 




w is\ Ke' 


J^d 

T m ~ J-foS 


s + 


Tm 


1 

rl + - 
*m 


'm 


s + 


*m 


(38) 


Error-network system: 




^is 
^ c + — 


1 + (T d + Ke r,S 


K e ' _ Td 
T d T m Jb s 


s 2 + 


1 JL 

Td + TjJ S + T d T m 


, 1 1 

+ | ■=; “ + -=- 
Td T m 


s 2 + 


T d T m 


+ 


*m 


T d T m, 


s + 


K e ’ 

T d T m 


Tachometer system: 


(39) 




w 


<Pr 


IS 


S + 


T t. 


'm 


Td 

JqjS 


+ 


r + i + fii 

T t T m T t T : 


s + 


T t T m 


K t ' 

T t ' T m T t T m 


s 3 + ( Jr + ir- + 


s 2 + 


1 Kg' A Kg' 
+ ) s + 


TtTm 


'm 


TtTm 


( 1 * 0 ) 
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An analysis using generalized response parameters can now be made 
if the closed- loop equations can be reduced by reasonable approximations 
to second-order equations having the form 








s 2 + 2£u n s + w n 2 


(4i) 


Equation ( 38 ) for the rate-gyro system is already in the desired 
form with T x = 0. If the assumption is made that Kg' » 1 (eq. ( 29 )), 
equation ( 38 ) is simplified to 


'Pb = 



V _?a 

T m " J b S 



s 


2 


+ 



s + 


K e ’ 

T m 



(42) 


Equation (39) for the error-network system has a second-order 
numerator and a third-order denominator; however, the numerator and denom- 
inator may be reduced by one order each by restricting the size of the 
lag time constant, t^, to small values. If it is first assumed that 


Ki* » 1 
T m » T d 


and 


T 


d 


« 



(43) 

(44) 

(45) 


and then the numerator and denominator of equation ( 39 ) are multiplied 
by Tq, the resulting equation is 


<Pb = 



t^s 3 + s 2 



s + 



(46) 


In the numerator the term Tqs 2 may be neglected with respect to s + l/r m 
up to frequencies somewhat above (j n , provided T d u n 2 « w n , or 


T d 



(47) 
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Similarly, in the denominator, the term T^s 3 may also he neglected, and 
equation ( ^-6) may he reduced to the desired form 


c Pb 



In this equation it is seen that 



and the assumption stated in equation (47) may he written as 


T d 



(48) 


(49) 


(50) 


The assumptions made in order to simplify the equation for the error- 
network system may he shown to he reasonable for most satellite control 
system applications. The assumption that Kj_' » 1 may he justified on 
the same basis that Kg’ » 1 is justified. The assumption that T m » 
is obviously valid since the motor time constant T m is large, and the 
inherent lag time constant is desired to he small. The validity of 

% « Kx’/Ke 1 may not he apparent from previous discussions; however, 
justification is readily provided by the combination of the derivative 
and straight-through paths for the error signal in sketch (c) for the 
error-rate network into a single lead- lag network with the following 
equation: 


Ki's 
1 + T d s 


1 + T d + 


Ki* 


1 + T d s 


(51) 


The lag time constant must he much smaller than the lead time constant 
and also much smaller than the reciprocal of the natural frequency 
(eq. ( 50 )) if the undesirable lag is to he insignificant during the 
transient response. 

Equation (4o) for the tachometer system also has a second-order 
numerator and a third-order denominator, which may he reduced by one order 
each by the same method used for the error-network system. If it is 
assumed that 


K e ' »lA t 


(52) 


and 


Kt’ » T m + T t 


(53) 
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and. if the numerator and denominator of equation (4o) are multiplied "by 
T t T m/ K tS the resulting equation is 


% = 




T-tKe* 

Kt* 


^d ( T t T m 2 
Jhs V^Kt' S 


+ s + 



T t T m 

K t ' 


+ s £ 


+ 


TtK e ' 

Kt’ 


s + 


K t ’ 


(54) 


In the numerator the term ( T t T m/Kt* ) s 2 and in the denominator the term 
( T t T m/- K t' ) s3 W neglected provided 


T t T m 

Kt’ 



(55) 


and equation (54) may he reduced to the desired form 


co 


+ 


is 
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Here, it is seen that 
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Wn J K t ' 

and the assumption stated in equation (55) may be written as 


T t T m 

Kt' 


« 


I W 

n/ K e ' 


(56) 


(57) 


(58) 


The assumptions made to obtain the simplified equation for the tachometer 
system are not quite as general as the ones made for the error-network 
system, and no attempt will he made to justify them. However, for any 
case in question, after values for the various parameters have heen 
selected, the simplifying assumptions may he easily tested to see whether 
the simplification was justified. 

Equations (42), (48), and (56) have the form of equation (4l) so it 
is a simple matter now to solve for the generalized linear parameters, 
w n , £ ) l/ T ik>n, and l/T 2 cj n , for the three systems. The results, including 
the assumptions used for simplifying the closed-loop equations, are listed 
in table I. The quantities l/^x^n and l/T 2 wn are used instead of 1 /t x 
and 1/t 2 because the response is affected by the size of l/id and 1/t 2 
relative to oj n . 
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The choice of the test numerical values for the generalized linear 
parameters is based mostly upon the following three considerations : 

(l) £ should he greater than 0.3 and less than 1.0 for good stability 
without sluggishness, (2) 1 /tx should be nearly equal to or greater than 
Wn so that the response for an initial displacement is not greatly 
affected by the position zero, l/Tj., which tends to make the response 
highly oscillatory, and (3) should be large enough for fast response. 

It should be mentioned here that the torque zero has a value 1/t 2 « w n 
for all practical cases. Hence, at the high frequencies of concern for 
transient response, where u > u n , the term s + i/t 2 reduces to s, and 
the numerator term associated with torque becomes T(p/Jb • Torque inputs 
are small and change only gradually, and therefore, they have little effect 
on the transient response. 


Numerical Example of Design Procedure 


The design criteria which have been established can now be applied 
to the numerical example of the Orbiting Astronomical Observatory (0A0) . 
From the most stringent fine control requirements presently contemplated 
for the 0A0, the following specifications for stellar sighting have been 
assumed: 

Maximum steady-state error: <Pessmax = second of arc 

Maximum initial error: cp eo = 60 seconds of arc 

Approximate settling time: ts = 1 minute 

Control time: t c = 105 minutes = 63OO seconds 

Maximum line-of-sight velocity: w ls = 0.01 second of arc/second 

Maximum roll (about line of sight) velocity: w r = 1.0 second of 

arc/ second 

Maximum disturbance torque: Tq = 100 dyne cm 

Body inertia: Jp, = 10 10 gram cm 2 

These specifications will now be applied to the selection of numerical 
values for all three systems. The outer loop gain, wheel angular momentum, 
and motor stall torque will be determined first. The outer loop gain must 
be large enough to satisfy equation (24) for w ls input and equation (27) 
for Ta_ input. For the 0A0 the gain requirement of equation (27) is much 
larger, and the resulting value for t c = 63OO seconds is 
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K e ' >130 (radians/sec) /radian 


(59) 


A value of K e ' with a slight margin of safety is now chosen as 
K e * = 150 (radians/ second) /radian. 

The required angular momentum capacity is readily obtained from 
equation (30) . When values for the OAO are substituted, 


Jv<jJw m ax = 6 . 3 ( 10) 5 gram cm 2 / sec 


(60) 


If the servomotor has a maximum speed of uv max = 6,000 rpm then the wheel 
inertia is Jv = 10 3 gram cm 2 . 

To determine motor stall torque (and the corresponding motor time 
constant) , equations (35) and (36) are plotted in figure 1. The solid 
line, representing equation (35) , shows that for a minimum settling time 
of 60 seconds, the motor stall torque must be at least 3 .3(10) 3 dyne cm. 

For this stall torque the motor time constant from equation (36) is 190 
seconds, which is large enough to produce nearly the same response as a 
T m of infinity. It should be emphasized that the motor stall torque 
given by equation (35) I s f° r minimum settling time with no overshoots. 

As will be verified later in the section on large- signal analysis, the 
actual stall torque required is nearly an order of magnitude larger than 
the minimum. For the present, a value of Tms is somewhat arbitrarily 
chosen as Tms = 2.1X10 4 dyne cm = 0.3 oz in. The corresponding value of 
motor time constant is T m - 30 seconds. This value of motor time con- 
stant is still large enough to have a relatively small effect upon the 
required motor stall torque as compared with the large effect of the 
overshoots in the response. 

Consideration will now be given to the linear design parameters, 

£, 1 /tj, and Wq. Using the first row of table I, it is seen that for the 
rate-gyro case, there is no zero at I/T3,. For the error-network and 
tachometer systems, it can be observed from table I that the following 
convenient relationship exists: 


t = 


T i w n 


(61) 


That several reasonable solutions exist for the error-network and tachom- 
eter systems is demonstrated by selecting values of £, such as 0.5 and 
0.7^ and calculating the corresponding values of l/T-jCJn, which are 1.0 
and 0.7, respectively. From figure 1.21 in reference 3 'the resulting per- 
centages of overshoot are found to be 30 percent for £ =0.5 and 20 per- 
cent for £ = 0.7. In the case of the rate-gyro system, where no position 
zero is present, the percentages of overshoot for £ = 0.5 and 0.7 are 
16 and 5 percent, respectively. 
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A damping ratio of £ = 0.7 will now be selected for all three 
systems, with the corresponding position zeros for the error-network an d 
tachometer systems located at I/t^ = 0.7. From formulas for damping 
ratio in table I and the values of T m and K e ' previously selected, the 
following are calculated: 

(62) 


and 

K 1 K * 1 / T 

or — i— = 2 t, / — ^ = 0.625 radian/ (radian/ sec) ( 63 ) 

rLe K -0 \ / -K-e 

When the above value of Kg' /he* is compared with the error requirement 
of equation ( 28 ) , it is found that 0.625 is satisfactory since it is 
considerably less than the maximum allowable value of 10. 

Since K e * has already been selected, Kg' and Kq' may be determined 
from equation ( 63 ), which gives Kg' or Kq* = 9k radians/radian. For the 
tachometer system, however, Tq and Kq' in equation (62) are not yet 
determined, since their selection is dependent upon the natural frequency, 
which will now be examined. 

For the rate-gyro and error-network systems, the natural frequency 
calculated from table I is w n = 2.24 radians/ second which seems reason- 
able for a settling time of one minute. If it were necessary, could 
be altered by changing either T m or Ke 1 , except that K e ’ must remain 
larger than the minimum given in equation ( 59 ). If a natural frequency 
of 2.24 radians/ second is used for the tachometer system, the resulting 
tachometer loop gain calculated from an expression in table I is 
Kt’ =30 radians/ (radian/ second) , and the value of T t calculated from 
equation ( 62 ) is Tq = 0.625 second. These values for Kq 1 and Tq must 
be held in abeyance until the simplifying assumptions are checked. 

Before the s im plifying assumptions for the three systems are checked, 
a value should be selected for Tq, which is the only remaining parameter 
to be chosen. The only demands on Tq are that it be small enough to 
satisfy the assumptions and prevent excessive delays, but large enough 
to be realizable. It would seem reasonable for a physical network that 
the lead time constant be at least ten times the size of the lag time 
constant. A value of Tq - 0.01 second will be selected, which is slightly 
over 60 times smaller than Kq'/Ke' . It will be shown later that the 
value of Tq can change over a wide range without significantly affecting 
the response. 

When all the values which have been chosen are substituted into the 
expressions in table I, the results shown in table II are obtained for the 
rate-gyro and error-network systems and Case A of the tachometer system. 


-=== = ■ = 0.114 sec 1 / 2 

jKt 7 n/K ^ 7 
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(Case B will "be discussed later.) A check of the simplifying assumptions 
shows that all of them are satisfied for the rate -gyro and error -network 
systems. For the tachometer system, however, equations (53) and ( 58 ) are 
not satisfied. Hence, equation ( 5 6 ) for the tachometer system is not 
valid when using the values chosen thus far. Further examination of the 
exact closed-loop equation is required to predict the response. Obtaining 
an exact solution is difficult because of the third-order denominator, and 
further analysis of this problem will be delayed until the large-signal 
response has been considered. 


LARGE- SIGNAL RESPONSE 


The linear analysis just completed is satisfactory for relatively 
small angular errors. However, for many satellites, the desired range of 
the attitude control system will be extremely large. In the particular 
case of the 0A0 the control system may have an initial error of up to 
60 seconds of arc, and still the system must control the vehicle attitude 
to within 0.1 second of arc. The criteria which have been used in the 
selection of K e ' imply that an error greater than 0.1 second of arc 
saturates the motor amplifier. Accordingly, a nonlinear analysis of 
system operation must be undertaken. 

Nonlinear Analysis 

The parameters which have been selected by the linear analysis for 
small-signal response, will now be used in a nonlinear analysis to obtain 
the settling time in a more exact fashion than the bang-bang study made 
earlier (eq. (35))- The large-signal settling time will be obtained, and 
methods for reducing this settling time to a minimum will be considered. 

A semigraphical approach called the " switching time" method, presented in 
reference 2, is well adapted to this problem. 

In order to apply the method, the block diagram for the three systems 
is first redrawn from sketches (b) , (c) , and (d) into the form shown in 
sketch (f) . The velocity and torque inputs are omitted from the sketch; 



Block diagram using loop 
gains and selector switch 


Sketch (f) 
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the loop gains are treated as single terms; the motor, wheel, and tody 
are combined into a type 1 * one, second-order plant; and a switch is included 
to select the method of damping. The plant is immediately preceded ty a 
box representing amplifier saturation where the level of saturation is 
given ty the fictitious voltage, which is defined as 


Vntaax _ KmVmmax “ “wftiax 


(64) 


For the systems shown in sketch (f), the following equations relating 
t sw to tAcp e and Acp e (where tAcp e is the time for the error to change from 
one point of zero velocity to the next, and Acp e is the incremental error 
change during tdcpg) are derived in reference 2 (eqs. (70) and (72)): 


Acp e - ax y^tgw - ■fcAcPeJ 

(65) 

/ \ 

■tAcpe = T m ln ( 2e Tm - l) 

(66) 


For these equations and the accompanying analysis to te valid, it is 
required that Vi%|rx « Ke'Acp e , so that the error signal saturates the 
amplifier most of the time. It is assumed that the duration of linear 
operation is negligible and that the reversal of full motor voltage occurs 
instantaneously as V' passes through zero. It is also assumed that the 
initial velocity (of the wheel or satellite) is zero; in many instances 
the initial velocity will te zero, and furthermore, analog computer results 
will te presented later to show that large initial velocities do not affect 
the general conclusions. 

The previous two equations can te solved simultaneously, and curves 
of t sw and tAcp versus Acp e can te plotted, using the value of "Vmmax 

for a particular system. For illustration the previous numerical example 
will te used. By substituting values which have already teen chosen into 
equation (64), it is found that Vjmax. = 6.3(10) “ 5 radians/second for all 

three systems. When this value for V m^ ay and % = 30 seconds are used 
in equations (65) and (66) , the solid-line curves of figure 2 are obtained. 


1 In this report "type” is used to refer to the number of integrations 

of the plant, or the number of poles at the origin (ref. 2 , p. 7) • 
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Response for minimum 
settling time 


Sketch (g) 



The usual application of these 
curves is illustrated in sketch (g) , 
which is very similar to sketch (e) . 

For cp e o - A<p e , the reversal of V' 
in sketch (f) at t sw causes the sys- 
tem to reach a condition of zero veloc- 
ity and zero position simultaneously 
at t^cp = ^Smin* Since there is no 
overshoot, there will he no further 
motion of the system output. For the 
large-signal analysis, it will he con- 
venient to use the curves of t S y and 
tAcp e In a slightly different fashion, 
as illustrated in sketch (h) . Here, 
cpeo ^ Acp e , and the error does not reach 
a condition of zero velocity and zero 
position simultaneously. However, 
since the motor voltage is reversed at 
the same time as in sketch (g) , the 
two curves still have the same shape. 
For known values of t sw and cp eo , the 
time for one-half cycle of operation 
is determined as tlcp e ; and the magni- 
tude of the overshoot is determined as 
Acp e minus cp eo . The overshoot in 
sketch (h) may he treated as the start- 
ing point for determining another half- 
cycle of operation. 


For the large- signal analysis, the 
rate-gyro system will he considered 
first. The relationship between t sw 
and cp eo ma 7 "be obtained by writing 
an equation for V' in sketch (f) . From the time of initial displacement 
until the first reversal time, the equation for V' for the rate-gyro 
system for a step input cp c ( s) = cp e0 /s (and the corresponding 
Vm'( s ) = vf iitoa.x/ s ) is 


Half cycle of response 
for 

Sketch (h) 


V’ 


'Pea 


Ke’ 


^max /K e ' + Kg' s \, 
s 2 \ 1 + T m s J 


(6?) 


s 
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Taking the inverse Laplace transform, dividing "by K e ' , and rearranging, 
we obtain 


V' 

K e ' 




t 


(68) 


Setting V' = 0 for t = tsw> and replacing (r m - Kg'/K e ' )/T m by the 
letter C, changes the equation into the approximate relation 


Teo 


= t V 1 

m v mmax 


j sv 

T m 



(69) 


This equation is used to determine a curve of tgv versus cp eo in 
figure 2, using a value for C = 0.979; which corresponds to the previ- 
ously selected values of T m = 30 seconds and Kg'/Ke' = O.625 
radian/ ( radian/ second) . If now an initial error angle of cp eo = 60 seconds 
of arc is assumed, then from the curve of tsw versus cp eo it is seen that 
t sw = 17-7 seconds. The corresponding value of Acp e from the curve of 
tsw "versus Acp e is 86.4 seconds of arc, and the value of t^cp from the 
curve of tAcp versus Acp e is 28.8 seconds. As illustrated in figure 3, 
the total angular change of 86.4 seconds of arc produces an overshoot of 
26.4 seconds of arc, and hence, the input for the next half cycle of the 
oscillation is 26.4 seconds of arc. Using cp eo = 26.4 seconds of arc in 
figure 2 gives t sw = 11.2 seconds, which results in tAcp e = 49* 2 seconds. 
The corresponding overshoot is found to "be 10.1 seconds of arc. The proc- 
ess can "be continued to as small an overshoot as desired, and the settling 
time is found "by adding the times for the successive half-cycles of oscil- 
lation. For the example "being considered, the settling time from cp e = 60 
seconds of arc to 0.1 second of arc is ahout 83 seconds as illustrated in 
figure 3. This value "becomes one point on a curve of t s versus cp eo in 
figure 2 (where outer scale is used for t s ) . Repeating this process for 
different initial errors gives values for the entire curve for t s . In 
figure 3 the motor drive voltage has been sketched to show the points at 
which switching occurs. 

It is interesting to note for comparison with sketch (e) that in 
sketch (g), which was obtained from figure 2, tsw is somewhat longer than 
one-half ts min » This actually indicates that a higher torque is applied 
during the deceleration period than during the acceleration period. The 
explanation for this is that the motor speed is increasing until tsw; as 
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a result , the reversal torque is larger than the stall torque as 
illustrated hy the "broken line in sketch (i) , whereas for sketch (e) the 
torque is independent of wheel speed. 


T m 



Sketch (i) 

The switching-time method can also "be used to determine an optimum 
nonlinear relationship "between gain and error signal which will reduce 
ts to "tAcp e = 'tsmin (with no overshoots) . The process will "be illus- 
trated for the rate-gyro system. The first step is to solve for K e ' in 
equation (68) with V' = 0 for t = t sw . 


K 


1 

e 



9 


eo 


T m w Wmax 


+ 1 - e 


Tsw 

Tm 


'SW 


‘m 


(70) 


The next step is to identify the error signal "by first writing the Laplace 
transform equation from sketch (f) for a step input. 
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_ ‘Peo 
Pe s 


^%ax 


(71) 


S 2 (l + T m s) 

Then taking the inverse transform and replacing t by tsw produces 

t< 

J YV » 

+ 1 - e “ - is I (72) 


<P< 


= T m V ' 


q> 


eo 


'SW 


m %iax tV 1 

' T m v i%iax 


‘m 


A comparison of the above equation with the denominator of equation (70) 


shows that 


\ *■ 
m\ 






1 - e 


J SW 


"SW \ 

T m 


(73) 


In the above derivation K e * was considered to be a constant; however, 
since the system is always limited except at the switching point, the 
value of K e ' is of no significance (as long as it is sufficient to cause 
saturation) except at tsw Hence, equation (73) may equally well be con- 
sidered as defining K e ’ in terms of instantaneous error, cp e , instead of 
cp e ]^ s w-’ since K e ' will then always have the desired value at t S w* In 

order to prevent overshoots the desired relationship between tsw and cp eQ 
may be obtained from the curve of t sw versus Acp e with cp eo = Ap e . 

For the numerical example being considered. Kg' = 94 radians/ radian 
and % = 30 seconds, giving a value for Kg'/T m of 

3*13 Cradians/ second) /radian. To illustrate how one point on the nonlinear 
curve of optimum gain versus error is obtained, a t sw of 13-5 seconds is 
arbitrarily chosen, and the corresponding input is read from the curve of 
tsw versus Acp e in figure 2 as cp eo = Acp e = 55 seconds of arc. Equa- 
tion (70) is then solved for Ke', which is found to be 
21.4 ( radians/ second) /radian, and equation (j2) is solved for cp e )t S w> 
which is found to be 20.7 seconds of arc. These values are used to plot 
one point on the curve of optimum gain in figure 4(a) (where cp e ]-t; has 
been replaced by cp e as previously justified) , which shows the relation- 
ship between outer loop gain and instantaneous position error required to 
produce minimum settling time for any initial error. 

The actual output-input relationship of the error amplifier producing 
K e ' is shown in figure 4(b) , where the ordinate is found by merely multi- 
plying K e ' by cp e from figure 4(a) . The constant gain prescribed by the 
previous linear analysis is illustrated in figure 4 by straight lines. It 
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should he mentioned that an optimum relation between Kg' and cp e to 
produce minimum settling time with constant K e ' could he calculated in 
the same manner. For the rate-gyro system the nonlinearity in K e ' 
required for optimum switching to produce minimum settling time may con- 
veniently he considered to represent a nonlinear sensor gain K s , since 
Kg' does not include K s and therefore the rate damping is not affected 
hy sensor gain. 


The error-network system will now he analyzed for large-signal 
response in the same manner as the rate-gyro system. The Laplace trans- 
form equation for V' from sketch (f) for a step input is 


V’ 


^eo 

s 


^max 

s 2 (l + T m s)_ 



+ 


Kj's \ 

i + r dS y 


( 7*0 


Taking the inverse transform, dividing hy 


K e ' , and rearranging gives 




Kj' 

T dV 


e 



T m v mmax 




TdKj' \ 

Tm KeV 



( 75 ) 


Letting V 1 = 0 for t = t sw , assuming T d « T m , and replacing 

( T m - Kd'/K e ')/' r m hy C changes the equation into the approximate 

relation 


Teo - T m Tmiax 


j sw 
T m 


- C 



( 76 ) 


which is identical to equation ( 69 ) for the rate gyro. Hence, the same 
curves may he applied, as shown in figures 2, 3> and except that Kq' 
is used instead of Kg' . It should he noted that unlike the rate-gyro 
system, the nonlinear gain in figure k cannot he applied to the sensor 
for the error-network system, since Kq' also includes K s and the rate 
damping is affected hy sensor gain. 
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For the tachometer system, the equation for V' is 


V’ 



Vm' 

m max 

s 2 



H" 


gtlgf 

1 + T t S ! 


1 + T m S 


(77) 


Taking the inverse transform and dividing hy K e * , 


y» 

‘Peo = T mV’n^ iax 





(78) 


Letting Y 1 = 0 for t = t sw . 


<Peo - 



(79) 


Switching- and settling-time curves for the tachometer system are 
shown in figure 5* For case A, Acp e is significantly larger than cp eo 
for small inputs. The resulting oscillatory response and long settling 
time are illustrated in figure 6, which was obtained hy the same graphical 
procedure described previously for the rate- gyro system. When the 
tachometer- loop gain, KtS I s increased by a factor of 100 and the damping 
ratio is kept nearly the same by increasing Tj. to 5 seconds, as shown 
in case B of table II, then the large-signal response becomes much more 
reasonable as shown in figure 6. It should also be observed that for 
case B in table II, the simplifying assumptions are satisfied. For the 
graphical computation of settling time in figure 5, tg was assumed to 
be approximately equal to t/^p for cp eo < 13 seconds of arc, because 

the method of reference 2 applies only for Acp e > cp eo . However, the 
approximation is reasonably valid in figure 5 because it is used for only 
a small portion of the settling time. In passing, it should be mentioned 
that the increased value of Kt* can no longer be completely neglected in 
equation (23) for the steady-state error due to Tq_. 

Obtaining the optimum gain requirements for reducing the settling 
time to a minimum for the tachometer system involves certain restrictions 
outlined in reference 2 and will not be attempted here. 
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Analog Computer Results 

In order to provide a check of the analyses previously presented, 
studies of wheel control systems were carried out on an analog computer. 
Using the same numerical example, small-signal simulation of the exact 
systems (third-order equations) verified the response predicted by the 
approximate linear analysis. In the presence of 100 dyne cm of external 
torque and a small initial displacement, all three systems tracked to 
within 0.1 second of arc for 105 minutes with a dynamic response and error 
buildup that agreed very well with that predicted by the second-order 
equations using generalized linear parameters. Only large- signal data 
will be discussed here because the linear response is well defined by 
generalized second-order curves. For the large-signal computer analysis 
the block diagrams of sketches (b) , (c) , and (d) were converted to the 
computer circuit shown in figure 7* 

The sensor output was limited at the values of cp Q , shown in 
figure 8. It is interesting to note that a limited sensor output results 
in a nonlinear sensor gain which is similar to the optimum gain shown in 
figure It-. Although the largest value of error limiting used was 
cp el = 10 seconds of arc, it was shown that very little change in response 
would have resulted from further increases in the size of the limit. 

Transient response curves for the three systems, using the parameters 
selected in the previous linear and nonlinear analyses, are presented in 
figures 9 and 10 . These curves represent large-signal response, and they 
correspond very well with the predictions made in figures 3 and 6. Graphs 
showing the effect of error limit on settling time are shown for the three 
systems in figures 11 , 12 , and 13 . From the figures it is apparent that 
settling time is essentially independent of the linear sensor range for 
all three systems for cp el > 6 seconds of arc. Some comparison can be 
made between figures 11 , 12 , and 13 from the analog computer and fig- 
ures 2 and 5 from the large-signal analysis. For the rate-gyro and error- 
network systems, the computer settling time of 85 seconds agrees quite 
well with the 83 seconds calculated settling time from figure 2 . For the 
tachometer system (ease B) the computer settling time of 75 seconds is 
somewhat larger than the calculated value of about 60 seconds from fig- 
ure 5, due to the approximation of t s by t^cp 0 for cpeo < 13 seconds of arc 
in the calculations. 

The dip in the settling-time curves for the rate -gyro ease of 
figure 11 apparently results from the curves of figure 8 approaching the 
optimum curves of figure 4 , thus the overshoot is reduced so that the 
settling time approaches -f&cp . It should be observed, however, that the 
dips in the curves of figure 11 still remain somewhat above the minimum 
value for settling time of tspj^ = lAcp e - 24 seconds shown in figure 2 
for &p e = 60 seconds of arc. For small error limits, the settling-time 
curve for the error-network system shown in figure 12 does not have a 
dip like that of the rate-gyro system, because the damping of the error - 
network system is reduced by limiting the error signal. For the tachometer 
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system there is a significant increase in settling time for small error 
limits as shown in figure 13 , which can "be explained from figure 5. For 
small effective values of cp eo produced by small cp ei , 9eo is greater 
than Acp e for a significant portion of the settling time, resulting in 
overdamped operation and long settling times. 

The effects on settling time of the time constants of the error- rate 
and tachometer networks are shown in figures 14 and 15. For the error- 
network system, t s is nearly independent of for Tq < 1.0 second. 

For the tachometer system, however, there is a broad minimum. When 
adjusting T-t to obtain this minimum, it must be remembered that the 
damping ratio for the linear region is directly proportional to the 
tachometer network time constant, t^. 

The graph in figure 1 6 of settling time versus motor time constant 
obtained from the computer results for the rate-gyro and error-network 
systems shows an approximately linear relationship. For the particular 
parameters selected, the settling time is about three times the motor time 
constant . 

The effect of initial wheel velocity on transient response is shown 
in figure 17 for the rate-gyro and error- network systems. Comparison of 
the curves shows that an initial wheel speed of one-half the maximum 
decreases the overshoot and settling time somewhat when applied in one 
direction, but increases the overshoot and settling time when applied in 
the other direction. Computer results were also obtained to show that 
initial body velocities of ± ( l/2) Wwmax^w/Jb = ±3»15 x 10~ 5 radians/second 

had similar effects on the shape of the response curve and lengthened the 
settling time to about 120 seconds. 


THEEE-DBffiNSIOE[AL PROBLEMS 
General Considerations 


Speculation about three-dimensional control points to at least two 
causes of unwanted cross-coupling torques. One of the most obvious causes 
is the change in direction of the angular momentum vector of a gyro or 
inertia wheel, which is well understood as a means of producing a torque 
about an output axis which is perpendicular to the input axis. Another 
apparent cause of cross-coupling torques is the misalinement of the wheel 
axes with the principal body axes, or the addition of an unbalancing mass 
between the axes of an otherwise balanced system, which results in pro- 
ducts of inertia being present in the body. Both of these causes were 
considered on the analog computer for the numerical example being con- 
sidered; however, before the results are described, some mention should 
be made of the angular momentum transfers required for satellite fine- 
control systems. 
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For a stellar-pointing satellite, such as the 0A0, any initial 
angular momentum of the tody is ah sorted into the wheels, and only small 
changes of wheel angular momentum are then required to overcome external 
disturbance effects. For an earth-pointing satellite any initial angular 
momentum of the body is absorbed into the wheels in the same manner as 
for a stellar-pointing satellite; however, large exchanges of angular 
momentum are required between the two wheels whose axes are in the orbital 
plane to maintain the total angular momentum constant as the satellite 
rotates in orbit. Thus the speed of two of the wheels will have an 
undamped oscillation at the orbital frequency. In general the control 
system response is relatively rapid, and the oscillation can be considered 
to be quasi-steady state. Three-dimensional design and simulation for an 
earth-pointing satellite, then, can be handled the same as for a stellar- 
pointing satellite by using a larger and slowly changing velocity input 
( u is) • 


Analog Computer Results 


Some analog computer results were obtained to study the theoretical 
three-dimensional behavior of the numerical example. The response was 
examined both with and without products of inertia for a system using three 
inertia wheels with mutually orthogonal axes. Tight control was maintained 
about the two axes which were perpendicular to the line of sight, and an 
undamped sinusoidal velocity was allowed about the line-of-sight axis with 
a maximum of u r = 50 seconds of arc/second, which is 50 times the value 
given in the specifications for the 0A0 example. 

The cross-coupling terms due to the directional changes of wheel and 
gyro momentum vectors were found to be relatively insignificant. Further- 
more, when inertia unbalance ranging from 1 to 5 percent of the body 
inertia was placed between the axes (to produce products of inertia) , its 
effect on the dynamic short-term response was negligible. When the roll 
rate of 50 seconds of arc/second was combined with near maximum wheel 
speed of several thousand rpm, the error exceeded the allowable maximum 
after several minutes. At the normal roll rate of 1 second of arc/second, 
however, the control time would not seriously be reduced below its 
estimated value of 105 minutes. 


EXPERIMENTAL TESTING 
Description of Apparatus 


To verify the calculations for design and to investigate the physical 
problems to be encountered, an experimental apparatus was assembled. Con- 
trol systems were mounted on a low- friction, three-degree-of-freedom 
platform, and control was maintain about the pitch and yaw axes, while the 
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platform was allowed to drift slowly about the roll axis. Data for the 
error-network system will be presented here, and some comment will be made 
about results obtained for a rate-gyro system which was also tested. The 
tachometer system was not tested because of the unavailability of any sort 
of tachometer device having the wide dynamic range demanded. A suitable 
device would have to be able to sense changes in speed of a fraction of 
an rpm while the wheel is turning at several thousand rpm. 

The main piece of test equipment for obtaining experimental data was 
a platform supported on an air bearing (fig. l8) . The actual air bearing, 
which is obscured in this picture by the platform, consists of a finely 
machined hall mounted in a cup with air injected from the bottom at about 
40 psig. The air bearing provides a low-friction support to simulate 
frictionless space conditions, and the 6-foot -diameter platform simulates 
the inertial mass of the body and provides an attachment surface for 
control- system instruments. Axl angular error detector was desired which 
would produce an output error signal essentially equivalent to that of a 
star tracker ; that is, a device whose output was sensitive to angular 
errors of fractions of a second of arc and linear over a region of many 
seconds of arc. The device selected is shown in sketch (j), including the 
light source and light sensor. Although this light source was uncollimated 
and has spectral characteristics, different from that of a star, it was 
adequate for this study. The output from the light source is obtained by 
using a ball bearing in a black box to reflect light which has been focused 
by a lens and filtered through a hole. Direct current lamp excitation is 
used to maintain constant light intensity. 



Light source and light sensor 


Sketch (j) 

The light from the light source is focused onto a wedge for one- 
dimensional control, which splits the light between two photomulti- 
pliers. An error signal is obtained by connecting the outputs of the 
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photomultipliers to a differential amplifer. For two-dimensional control, 
the wedge is replaced hy a four-sided pyramid and two additional photo- 
multipliers are used. The region of linear sensor operation is determined 
"by the spot size on the pyramid, since the voltage output is constant when 
all of the spot is focused onto one pyramid face. 

The d-c error output from the sensor (i.e., from the differential 
amplifier) was modulated "before it was amplified "by a-c amplifiers, which 
were used to minimize drift. Since a-c lead networks posed operational 
problems, the lead signal for the error-network system was obtained by 
demodulating, d-c shaping, and remodulating. The motors used to drive the 
inertia wheels were a-c servomotors. 


Inputs and Parameters 


After pendulosity torques were minimized by shifting vertical weights 
and the platform was finely balanced (within 1,000 dyne cm) by shifting 
horizontal weights, the disturbance torques were measured. The three most 
important torque inputs were the air-bearing turbine torque, the air-mass- 
circulation torque, and the building-vibration torque. Magnetic torques, 
which were already very small, were minimized by constructing the platform 
from nonmagnetic materials. The turbine torque, which acted about the yaw 
axis, was due to imperfect air flow in the bearing seat, and seemed to be 
relatively constant at nearly 10,000 dyne cm in the counterclockwise direc- 
tion. Random air mass circulation and building vibration torques produced 
peak inputs on the order of 10,000 dyne cm, acting predominantly about the 
pitch and roll axes. Position error inputs up to roughly 60 seconds of arc 
were obtained by blowing gently on the platform. 

It should be mentioned that position measurements which were taken 
with the sensor previously described 2 were accurate only in a dynamic 
sense, since the sensor was used to measure accurately position relative 
to its null signal position, and no attempt was made to minimize slow 
drifts in the sensor in order to measure absolute pointing position 
relative to the actual direction of the light source. 

The parameters used for the error-network system may now be presented. 
The platform inertia had a value of Jq = 3XIO 8 gram cm 2 , and the wheel 
inertia had a value of J w = 7*8xi0 3 gram cm 2 . The servomotor and inertia 
wheel combination had the following characteristics: This = 0.6 in. oz, 

Vnijjiax = 1*0 volts, T m = 100 seconds, and Ww max = 560 radians/second. The 

resulting motor gain was K m = ^Wmax/^mmax = ^ ( radians/ second) /volt . 

Since the external torque inputs were high, the motor amplifier was made 
to saturate for an error of 10 seconds of arc instead of 0.1 second of arc. 
The specific gains used were Ks = 1.2X10 2 volts/radian, K e = 6.7 
volts/volt, and Ka = 100 volts/volt. Using these values the outer loop 

2. 

These measurements were also verified with a precision autocollimator. 
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gain was found to be K e ' = 290 ( radians/ second) /radian. A damping ratio 
of £ = 0.7 was selected^ which required that the error-network loop gain 
have a value of Kq 1 = 240 radians/radian, or Kq = 55O volts/ (volt/ second) - 

The parameters which were selected can be used to predict the small- 
signal and large-signal responses. As mentioned previously; a damping 
ratio of £ = 0.7 results in a position zero at I/t;^ = 0.7 and a cor- 
responding overshoot of 20 percent. From the expressions in table I, the 
natural frequency is found to be w n = 1.7 radians/second; and the zero for 
torque input is found to be l/i^wn = 0.0061. A response curve obtained 
on the analog computer for a small-signal initial error of cp eo = 10 sec- 
onds of arc is shown in figure 19(a) . The indicated overshoot of 20 
percent agrees with the predicted value. 

The large-signal response can be roughly predicted by the switching 
time method. Values for tsw versus Acp e and tAcp e versus Acp e obtained 
from equations (65) and (66) are plotted in figure 20. The curve of 
tsw versus cp eo comes from equation (76) ; where the value of C is found 
to be 0.992. It should be observed in figure 20 that for cp eo < 78 sec- 
onds of arc; Acp e is less than cp eo . Hence; for a completely nonlinear 
system; no overshoot would be obtained and the response would be sluggish. 

It can be observed that the minimum possible settling time for cp eo = 60 
seconds of arc is about 3 seconds. The analog computer results presented 
in figure 19(b) show that the actual settling time is slightly over 
5 seconds . 


Experimental Results 


The actual experimental results obtained for the platform response 
are presented in figure 21; which compares quite well with the analog com- 
puter results of figure 19. As shown in figure 21; transient response 
during the settling period was fast with very little overshoot; and tight 
control was maintained after the settling period. A long time history 
showed that the random error remained less than ±0.5 second of arc. This 
random error was probably caused by the random torque inputs. The small 
ripple in figure 20(a) is due to electrical noise. The larger overshoot 
and longer settling time in figure 2l(b) compared to figure 19(1) is 
probably due to initial wheel velocity. The other minor dissimilarities 
between figures 19 and 21 were not accounted for; however, some general 
comments about practical problems which cause differences between theory 
and experiment may be made . 

Any unbalance or pendulosity remaining in the platform after balancing 
produces a steady torque which may aid, but usually hinders, tracking. 
Misalinement of the sensor relative to the center of platform rotation 
causes cross coupling, so that the roll motion allowed about the line-of- 
sight axis couples through the other two axes and changes the pointing 
direction of the line-of-sight axis. Electrical delays add small lag 
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terms to the open- loop equations , hut they are generally negligible for 
the systems considered here. The effects of saturation in the sensor and 
in the motor amplifier have already been discussed. Although hysteresis 
and dead zone have no obvious effects in the experimental results for the 
error-network system, they deserve some consideration. Some computer 
results were obtained by modifying the circuit of figure 7 to include 
hysteresis and dead zone in the sensor, gyro, and motor amplifier. It 
was found that hysteresis in any of these components caused limit-cycle 
response. Furthermore, limit cycles resulted from dead zone in the gyro. 
A dead zone in the sensor helped to stabilize the system, but it resulted 
in a steady-state error as large as the dead zone. A dead zone in the 
motor amplifier had very little effect on the transient response, but it 
also produced a steady-state error as large as the dead zone . 

The results for the rate-gyro system should also be mentioned. The 
large-signal response was essentially the same as for the error-network 
system, as predicted in the analytical discussion. However, the small- 
signal response was dominated by limit-cycle operation, which resulted 
from the rate-gyro dead zone. The threshold of the gyros was about 
30 seconds of arc/second, and the error amplitude of the resulting limit 
cycle was ±5 seconds of arc. 


CONCLUSIONS 


To design some reasonable systems for precise attitude control of a 
satellite, several methods of damping were first considered to provide the 
necessary compensation. It has been shown theoretically that either a rate 
gyro or an error-rate network is suitable without modification, but that a 
tachometer must be connected in series with a high-pass filter in order to 
keep the error signal from building up with wheel speed. 

The effects of adding an integrator to the control loop in order to 
minimize the error have been considered, and it has been shown that, for 
the requirements of the Orbiting Astronomical Observatory and other similar 
applications, the use of an integrating gyro (or electronic integrator) 
requires a much larger motor than would otherwise be necessary. 

Of the various possible control systems considered, three were 
selected for a detailed analysis. The three systems, which differ only by 
the method of damping, are referred to as the rate-gyro, error-network, 
and tachometer systems; the tachometer system includes a high-pass filter. 

An analytical investigation, which consisted of a generalized linear 
analysis, a nonlinear analysis using the switching-time method, and an 
analog computer study, showed that all three systems were theoretically 
capable of producing adequate response and also of maintaining the required 
pointing accuracy for the Orbiting Astronomical Observatory of ±0.1 second 
of arc. Practical considerations and an experimental investigation using 
laboratory apparatus showed, however, that the error-network system was 
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superior to the other two. The rate-gyro system was shown to he inferior 
because the threshold level causes a significant amount of limit-cycle 
operation. The tachometer system was shown to he inferior because a device 
with the required dynamic range of operation does not appear to he avail- 
able; the response is sluggish for low limits on the sensor output; and 
the large amount of tachometer feedback had a significant effect on the 
pointing error due to disturbance torque. For the error-network system, 
experimental data taken under laboratory conditions with relatively large 
extraneous disturbances show that a dynamic tracking error of less than 
±0«5 second of arc was obtained. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Nov. l6, i960 
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APPENDIX A 

VELOCITY AND TORQUE INPUTS 


In this appendix the velocity and torque inputs that will he 
encountered during stellar observations by the Orbiting Astronomical 
Observatory are briefly described. 


VELOCITY INPUTS 


The line-of- sight velocity inputs consist of parallax and aberration 
components which result from the motion of the satellite relative to the 
earth and from the motion of the earth relative to the sun. 


Velocity aberration occurs because the velocity of the satellite is 
a significant fraction of the velocity of light, as illustrated in the 
sketch. The aberration angle is a maximum when the satellite velocity is 



Velocity aberration resulting from orbital motion 
of a satellite about the Earth 


Sketch (k) 

perpendicular to the light velocity. For a 500-mile circular orbit having 
a period of 100 minutes, the satellite velocity is 4-72 miles/second. The 
maximum aberration angle found by taking the ratio of satellite velocity 
to the velocity of light is 2.54X10" 5 radian or about 5 seconds of arc. 
Since the aberration angle is sinusoidal at orbital frequency, the maximum 
line-of-sight velocity is found by multiplying the maximum aberration angle 
by the angular frequency, which gives 1.6X10 -6 radian/minute = 0.005 second 
of arc/second. This component of line-of-sight velocity due to aberration 
is very much larger than the component resulting from the motion of the 
earth about the sun, which may be completely neglected. 


Parallax, which is the difference in direction of a body when viewed 
from different places, may also be completely neglected, because of the 
large distances which separate the earth from the stars. 
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TORQUE INPUTS 


The disturbance torques acting on an earth satellite result largely 
from the earth's gravity gradient, the sun's radiation pressure, the 
earth' s atmospheric density, and the earth' s magnetic field. Only the 
order of magnitude of these torques will be discussed here; a more detailed 
discussion may be found in reference 1. The torque due to the earth's 
gravity gradient has been frequently discussed. For the astronomical 
satellite, a one percent inertia unbalance could result in a torque whose 
amplitude varies at twice orbital frequency between zero and 160 dyne cm 
with an average value of 80 dyne cm about one axis. 

Solar radiation causes a significant torque if the body is 
geometrically unsymmetrical or if the surface materials have different 
reflectivities so that the pressure center is not in line with the mass 
center when viewed from different angles. Using a solar pressure of 
7X10" 5 dyne/cm 2 for a body of medium reflectivity, a torque of 100 dyne cm 
will be produced on an exposed surface area of 100 square feet if the 
pressure and mass centers differ by about 15 cm or l/2 foot. 

Atmospheric density torque may be treated in the same manner as solar 
radiation, since both effects are calculated using the distance between 
the center of pressure and the center of mass. For the 500-mile orbit 
being considered, the dynamic pressure is roughly 1.5X10" 5 dyne/cm 2 , 
resulting in a torque which is somewhat smaller than the solar torque. 
However, if the altitude were reduced to ^00 miles, the atmospheric torque 
would increase by a factor of about seven, so that it might well be the 
largest disturbance torque. 

When the earth' s magnetic field interacts with the net magnetic field 
of the satellite, another unwanted torque is produced. Since the field of 
the satellite is produced by all of the magnetic materials and circulating 
loop currents distributed throughout the satellite, the resulting torque 
is very difficult to estimate. Although large torques could result from 
this source, they may be significantly reduced by nonmagnetic design and 
symmetrical distribution of components. 
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APPENDIX B 

ANALYSIS OF AN INTEGRATING -GYRO SYSTEM 


In this appendix an integrating -gyro system will he compared to the 
rate-gyro, error-network, and tachometer systems. It will he shown that 
although a lower outer loop gain is allowable, a much larger motor is 
required hy stability and error requirements for the Orbiting Astronomical 
Observatory and other similar applications. For this analysis the closed- 
loop response equation and the closed- loop error equation will be examined. 
These equations for the integrating gyro may be obtained directly from 
those for a rate gyro by replacing K a by K a /s. With this modification, 
equation ( 38 ) for closed- loop response becomes 


Tb = 



(Bl) 


Applying Routh' s criterion shows that the system is stable if 

T <^- 

Tm < K e ' 


(B2) 


In order to develop error criteria, the closed- loop error equation 
for an integrating gyro is obtained by replacing Ka by K a / s in equa- 
tion (4), giving 


9e = 




s 


3 




s + 



(B3) 


The steady-state errors for step inputs are found by applying the final 
value theorem. 

For cp c (s) = cp c /s (with cj 1s = Tq = 0) , 


^ess 

<P C 


lim scfu = 0 
s-O 


(b4) 
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For w is (s) = w is /s (with cp c = Tq = 0) , 

9ess _ Kg' 
“is ~ K e ' 

For Td( s) = T s (with cp c = m is = 0) , 

9ess _ _1 1 

Td Jb V 


(B5) 


(b6) 


To establish one error criterion , 


equation (B5) can be written as 



Tess- 


'max 


w is 


(B7) 


Incidently, this is identical to equation (28) for the rate gyro. For the 
Orbiting Astronomical Observatory the maximum ratio resulting from 
cp esSmax = 0.1 second of arc and ojis = 0.01 second of arc/second is 

Kg' /K e ' = 10 radians/ (radian/ second) . 

To establish another error criterion, equation (B6) may be written 
as 


C ' > — 

6 - Jb 9 


ess max 


(B8) 


The minimum allowable outer loop gain resulting from 9 ess max = secon< 3- 
of arc = 5X10 -7 radian, T^ = 100 dyne cm, and Jq, = IO 10 gram cm 2 is 
K e ' = 0.02 (radian/ second) /radian, which is much smaller than the outer 
loop gain required for the rate-gyro, error-rate network, or tachometer 
systems. 

The stability requirement of equation (B2) may be combined with the 
error requirement of equation (B7) to determine the maximum allowable 
motor time constant. If Kg'/Ke' is chosen as its maximum of 
10 radians/ (radian/ second) to allow T m to be as large as possible, then 
by equation (B2) , T m must be less than 10 seconds. When the angular 
momentum due to disturbance torque has been determined from equation (30) , 
the motor stall torque may be calculated by using equation (38). If J-w 
is selected as before as 10 3 gram cm 2 for uv max = 6000 rpm and T m is 

less than 10 seconds, then the stall torque may be written as the inequal- 
ity .> Tims > 6.3XIO 4 dyne cm, or T ms >0.9 oz in. It should be emphasized 
that the minimum value of stall torque indicated here corresponds to a 
neutrally stable system having sustained oscillations. Actually, analog 
results and a detailed analysis indicate that to have reasonable damping 
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ratios, the stall torque must range from 3 to 10 oz in. and is therefore 
10 to 30 times larger than the value of 0.3 oz in. determined for the 
rate-gyro, error- network, and tachometer systems. If weight and power 
must he conserved, the integrating-gyro system is not as promising as the 
others and therefore was not considered in detail in this report. 

It should he observed that for the numerical example just considered, 
the motor size for the rate-gyro, error- network, and tachometer systems 
was determined from the settling-time requirement while the motor size for 
the integrating-gyro system was determined from the steady-state error and 
stability requirements. It must be conceded that the integrating-gyro 
system might compare favorably with the other three systems if a much 
shorter settling time were required, which is not likely for precise sat- 
ellite attitude control, or if a much larger ratio of cp esSmax to w is in 
equation (B7) were allowed, which is entirely possible for some 
applications. 
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APPENDIX! C 

COMPLETE MOTOR AMD WHEEL BLOCK DIAGRAM 


A complete Block diagram of the motor-wheel combination for 
controlling satellite attitude is shown in the sketch. Here B e is the 



Sketch (l) 

hack emf constant, KiJ3 e i- s the forward voltage constant. Bp is a viscous 
friction constant, and Bjj is a magnetic torque constant resulting from 
the earth's magnetic field. The externally applied wheel torque, Tq w , is 
composed of such effects as gravity torques on an unbalanced wheel, and it 
is completely negligible. 

The wheel velocity, resulting from T m - Tq w is relative to 
inertial space, as also is the body velocity, resulting from T m - Tq. 
The wheel velocity relative to the body, Wy-q,, is actually equal to uv 
minus uq. However, since the body inertia is so much larger than the wheel 
inertia, w^b can be approximated by (%. 

Neglecting Tq w and approximating Wy-q by (%, the motor-wheel 
transfer function becomes, 


NmB e (B in + Jys) 
jEm _ Be + Bf + Bm 
V m J w s 


(Cl) 


B e + Bf + Bm 
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After the magnetic torque constant, B m , which is very small, and the 
viscous friction constant, which is relatively insignificant, are 
neglected the transfer function reduces to the following form: 

Tm _ KfflJ-yS 

Vm , J w 
1 + B e S 

This form is used in sketches (t>) , (c) , and (d) with J w /B e replaced t>y 
the motor time constant, T m - 

A useful equation for the motor time constant may he derived hy using 
the following approximate equation taken from the sketch: 


(C2) 


Bm = KmBeVm - Bguv 


(C3) 


By setting ov = 0 and V m = the stall torque is obtained. 


Bms - KjnBeVmmax 


(c4) 


and by setting T m = 0, V m = , and - GV max , equation (C3) gives 


Be^Wmax = K^eV^. 


.ax 


(C5) 


By combining equations (c4) and (C5) it is found that 


B e = 


Ti 


•ms 


^Wi 


(c6) 


max 


When this relation is used, the equation for T m becomes 

Jw^wmax 


T_ = 

m T , 


ms 


(C7) 


This form is used in equation ( 36 ) of the report. 
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TABLE II.- TABLE OF SUBSTITUTIONS FOR THE GENERALIZED LINEAR PARAMETERS AND 
SIMPLIFYING ASSUMPTIONS USING THE NUMERICAL EXAMPLE 


, j r\ t 

— or — L 


SYSTEM tu n (radians/sec) C 


2.24 


Er 

network 


Tachometer 
(Case B) 


2.24 



2.24 

.7 

.224 
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Simpliflying 
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K e K e ‘ Ke 
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\ radian/sec / T d or r t (sec ) 


.015 

94 » 1 
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Figure 1.- Minimum settling time and. motor time constant versus motor 

stall torque. 





Figure 2 


Incremental position error, A<£ e , and initial position error, <£ eo , seconds of arc 

- Switching- and settling-time curves for rate-gyro and error-network systems 
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Figure 3*“ Transient response and motor voltage curves for rate-gyro and error-network systems. 
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Figure 5*“ Switching- and settling-time curves for tachometer system 
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Figure 6.- Transient response 
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3 - T d /J b 


Scales: 


4 - K e /3.I25 

5 * K 0 K m /2xlO-" 

6 ~ IO/r m 

7 " 0.04 

8 - 2 x I0 6 J w /J b 

9 " 0.1 
10- Kg/2.5 

I I - Kj/6.25 

12 ~ O.OI/T d 

13 - 4xl0 5 Kt 


^eo 3x I0“ 6 radians /volt 

fie 3 x I0~ 6 radians/volt 

^ei 6 xIO' 7 radians/volt 

V 3,75 x I0" 7 volts/volt 

V m K m 75 radians/sec/volt 

u) b 1.5 x I0" 6 radians/sec/volt 

u) )s 1.5 x I0~ 6 radians/sec/volt 

T d /J b 7.5 x I0" 8 radians/sec 2 /volf 


l4"0.5/r t oj w 15 radians/sec/volt 

15 ‘ w w0 a» w0 15 radians/sec/volt 


Note: Double time scale was used 


Figure 7*- Analog computer circuit with selector switch for rate-gyro, error-network, and tachometer 

systems. 
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Figure 10.- Transient response curve for tachometer system as simulated ty analog computer 
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Figure 12.- Settling time versus linear range of sensor for error-netvork 

system. 
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Figure 15.- Settling time versus tachometer network time constant for 

tachometer system. 
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Figure 1 6 .- Settling time versus motor time constant for rate-gyro and 

error-network systems. 
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Figure 17.- Transient response curves shoving effect of initial wheel speed for rate-gyro and error- 

network systems as simulated hy analog computer. 
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Figure 18.- Laboratory apparatus for experimental testing. 
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(b) Large-signal response. 

Figure 19.- Small- and large-signal transient response curves for error 
network system as simulated "by analog computer . 




Switching time, t sw , and incremental error time, t , seconds 


3.6 


-4 

o 



Figure 20.- Switching- and settling-time curves for rate-gyro and error-network systems 
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(a) Small-signal response. 
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(b) Large-signal response. 

Figure 21.- Small- and large-signal transient response curves for error 
net-work system obtained by experimental testing. 
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